On the existence of a finite-temperature transition in the two-dimensional gauge glass 
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Results from Monte Carlo simulations of the two-dimensional gauge glass supporting a zero- 
temperature transition are presented. A finite-size scaling analysis of the correlation length shows 
that the system does not exhibit spin-glass order at finite temperatures. These results are compared 
to earlier claims of a finite-temperature transition. 
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The existence of a finite-temperature spin-glass transi- 
tion in the two-dimensional gauge glass has been a source 
of great controversy. While the existence of a transi- 
tion in three dimensions is undisputedji several differ- 
ent approaches&2iaiai2iMi2»12i have been made to deter- 
mine the spin-glass ordering temperature T c in two di- 
mensions with different results. As the gauge glass in 
two dimensions is a model commonly used to describe 
the vortex-glass transition in dirty high-temperature su- 
perconductor thin films, 11 it is important to settle the 
question of whether the system exhibits spin-glass or- 
der at finite temperatures or not. Experimental re- 
sults on YBCO thin filmsi^ show no vortex-glass or- 
der, whereas bulk materials do exhibit a transition to 
a glassy phase. A brief summary of the different ap- 
proaches used to determine T c is presented, followed by 
results from a novel technique^ that strongly supports 
the absence of a glassy phase at finite temperatures in 
the two-dimensional gauge glass. 

Early work by Fisher et al? who study the size- 
dependence of domain-wall energies with a zero- 
temperature domain-wall renormalization group method, 
shows evidence of a zero-temperature transition in the 
two-dimensional gauge glass. Gingrasji^ as well as Akino 
and Kostcrlitz^ who also study the stiffness exponent 
of the model, find agreeing results. 

In the aforementioned work, 2 Fisher et al. also ana- 
lyze the scaling of Binder ratios^ of the order parameter 
and obtain further evidence of a zero-temperature tran- 
sition. These results are in agreement with Monte Carlo 
simulations by Reger and Young£ who conclude that the 
transition temperature is not finite by studying the scal- 
ing of the currents induced by a twist in the bound- 
ary conditions. Katzgraber and Young 9 find evidence 
of a zero-temperature transition for the two-dimensional 
gauge-glass Hamiltonian with periodic boundary condi- 
tions by studying the scaling of the currents and the 
spin-glass susceptibility. They also conclude that one 
can scale data for the susceptibility well for any posi- 
tive T c , probably due to the extra fitting parameter rj, 
the anomalous dimension of the spin-glass order parame- 
ter q, and show that studying the scaling of the spin-glass 
susceptibility alone is not enough to obtain evidence of 
a finite-temperature transition. This may explain why 
similar simulations by Choi and Park^ who study the 
susceptibility only, find a finite-temperature transition 



with T c = 0.22 ± 0.02. 

Furthermore, Kim 8 has studied the gauge glass with 
fluctuating twist boundary conditions^ by means of 
resistively shunted junction (RSJ) dynamics and finds 
a finite T c in agreement with the work of Choi and 
Parkji and with early results by Lii& (challenged later 
by Simkm 19 ). Surprisingly, Granato,- who also uses RSJ 
dynamics but with periodic boundary conditions, obtains 
T c = 0. 

More recently Holme et al}° have studied the two- 
dimensional gauge glass by means of Monte Carlo simula- 
tions with fluctuating twist boundary conditions. By in- 
troducing an extra anomalous exponent b they are able to 
scale data for the currents and the helicity modulus and 
so obtain a finite transition temperature close to 0.20. 

As there are diverging opinions on the spin-glass tran- 
sition temperature of the two-dimensional gauge glass, 
the issue is reconsidered in the present work using a dif- 
ferent approach. The scaling of the finite-system correla- 
tion lengths is studied, giving compelling evidence that 
T c ss 0, at least for the system sizes studied here. These 
results are supported by data for the spin-glass suscepti- 
bility. 

The gauge-glass HamiltonianjSi which describes a dis- 
ordered granular superconductor, is given by 
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where the sum ranges over nearest neighbors on a square 
lattice of size N = L 2 and cf>i represent the angles of 
the XY spins. The Aij are quenched random variables 
uniformly distributed between [0, 27r] with the constraint 
that = —Aji. In this work J = 1 and periodic bound- 
ary conditions are applied. 

Previously^ we estimated the critical temperature of 
the two-dimensional gauge glass by studying the scaling 
of the currents^ induced by a twist in the boundary con- 
ditions. This has an advantage over studying the crossing 
of Binder ratios— because the latter will not splay enough 
for XY systems^ at T < T c , making it almost impossible 
to accurately determine the crossing point. Studying the 
currents requires some care as well, as these have differ- 
ent finite-size scaling forms^ depending on a finite or zero 
transition temperature. Finally, the currents show strong 
statistical fluctuations and so require a large amount of 
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FIG. 1: Correlation length £l divided by L for different sys- 
tem sizes. At all temperatures the data decrease with in- 
creasing L. This indicates that if T c is finite, it must be less 
than 0.13, the lowest temperature simulated. The shaded re- 
gion r epre sents the area around T ~ 0.20 where the claimed 
(Refs. l7lgUlCl) finite-temperature transition takes place. The 
inset shows a scaling plot of the data for / L plotted against 
L 1,V T for T c = and l/v = 0.39. One can see that the data 
scale very well, especially for low temperatures. 



disorder realizations. Therefore in this work the scaling 
of the correlation lengths is studied as these show little 
statistical fluctuations, have a simple finite-size scaling 
form and splay out well for T < T c . This method has 
been applied with success to Ising^ as well as vector spin 
glasses^ The spin-glass order parameter for the gauge 
glass is given by 
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where "a" and "/?" represent two replicas of the system 
with the same disorder. Equation J5J is generalized to 
finite wave vectors k: 
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The wave vector dependent spin-glass susceptibility 
Xsg ( k ) is then g iven by 



XsG (k) = 7V[(| g (k)| 2 )] E 



(4) 



Here, (• • •) denotes a thermal and [• ■ -] av a disorder av- 
erage. The finite-system correlation length £l is deter- 
mined from an expansioni&SiSi of the spin-glass suscep- 
tibility for small k£i in the framework of an Ornstein- 
Zernicke approximation 25 



Keeping only the leading term in the expansion and tak- 
ing into account the lattice periodicity as well as correc- 
tions to scaling, we obtain 
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where k m ; n = (2tt/L,0, 0) is the smallest nonzero wave 
vector. In Eq. @, x SG (0) is the standard spin-glass sus- 
ceptibility. 

The finite-size scaling of £l can be understood by 
studying the finite-size scaling of the spin-glass suscepti- 
bility extended to finite wave vectors k 

XsG (T, L, k) = L^>C[L^ V {T - T c ),kL] . (7) 

From Eq. © /c£ in £f, ~ X SG (0)/x SG (fe min ) up to an addi- 
tive constant, and so the L-dependent prefactor in Eq. Q 
cancels. The dependence of x SG on kL drops out as well, 
as in Eq. k = fc m ; n = 2tt/L. We therefore obtain the 
simple finite-size scaling form 



^ L /L = X[L 1 ' V {T-T C )] , 



(8) 



[XsgW/XscW]- 1 = 1 + (k£ L ) 2 +0[(k£ L ) 4 ] 



(5) 



and so, the introduction of an anomalous exponent^ is 
not plausible. This has also the advantage that the finite- 
size scaling has few parameters. From Eq. JBJ we see 
that data for different system sizes L will cross at T = 
T c . If T c — 0, one expects that data for different L will 
decrease with increasing L and not cross for any finite 
temperature. 

For the simulations the parallel tempering Monte Carlo 
method2S*21 is used. Because the equilibration test for 
short-range spin glasses introduced by Katzgraber et 
a £ | 22 J 28 4 29 (j oes no t WO rk for the gauge glassS as the dis- 
order is not Gaussian, equilibration is tested by the tra- 
ditional technique of requiring that different observables 
are independent of the number of Monte Carlo steps. 
By doubling the number of Monte Carlo sweeps between 
each measurement until the last three agree within error 
bars, equilibration is ensured. The equilibration times 
used for the different system sizes are listed in Ref. |jj 
Table I. For the largest size studied, £ max = 16, 10 3 dis- 
order realizations are performed, whereas for the smaller 
sizes L, 10 3 (L max /L) realizations. The lowest temper- 
ature studied is T — 0.13, well below the claimec&2ii£ 
transition at T w 0.20. Here, 30 temperature replicas 
are used with the highest temperature being T — 1.058. 

Figure^shows data for £l/L for different system sizes 
L. The data decrease with increasing L for all temper- 
atures. In particular, there is no crossing of the data 
around T w 0.20 (shaded area) where the spin-glass tran- 
sition is supposed to take placeiSiiS The inset of Fig. 
shows a finite-size scaling plot of the data according to 
Eq. |JSJ. One can see that the data scale very well for 
T c = and l/v = 0.39 ± 0.02, values found in Ref. H 
by scaling the currents and the spin-glass susceptibility. 
The error bar in l/v is determined by varying its value 
until the data do not scale well. An attempt to scale the 
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FIG. 2: Data for ln(x SG /£ 2 ) vs ln(L) for different temper- 
atures. At T c one expects Xsg ~ L 2 ~ v . The data show a 
strong downward curvature indicating that T > T c for the 
range of temperatures shown. This is emphasized in the inset 
where the curvature (from a second-order polynomial fit) is 
displayed. At T c one expects the curvature to cross zero. Here 
the data extrapolate to T c ~ (the dashed line represents a 
second-order polynomial fit extrapolation). The shaded area 
shows the region where one would expect the curvature to 
cross zero if T c ~ 0.2. A finite-size scaling analysis of the 
spin-glass susceptibility according to Eq. Q works well for 
low enough temperatures with T c — 0, 77 = 0, and l/u = 0.39, 
in agreement with results from Ref. 0- 



data in Fig. ^to any value of a finite T c , in particular 
0.20, fails for any choice of v. 

Standard finite-size scaling predicts that the spin-glass 
susceptibility scales a a 30 i 31 



X SG =L 2 ^C[L^(T~T C )} 



(9) 



which means that at T = T c it should scale as a power law 
- L 2 -v. In Fig.IUdata for \n( XsG /L 2 ) vs ln(£) is shown. 
If the data were at T c , one would expect a straight line. In 
particular, T = 0.20 shows a strong downward curvature 
indicating that the system is above the critical transition 
temperature. A fit of the data in Fig. [3 to a straight 
line for the different temperatures yields quality of fit 



probabilities 3 ^ between 10 -9 for the lowest, and 10~ 20 
for the highest temperature, respectively. A second-order 
polynomial fit to the data, instead, has quality of fit prob- 
abilities ~ 0.9. The curvature of the second-order fits is 
determined and shown in the inset of Fig. [5] Because 
at T < T c [T > T c ] one expects Xsg(^) to bend upward 
[downward], the curvature has to go through zero at T c . 
The data in the inset of Fig.|21do not cross the horizontal 
axis for the temperature range studied and extrapolate 
within error bars to T c « 0. This is in agreement with 
results by Katzgraber and Young 9 ^ who obtain good scal- 
ing of the susceptibility at low temperatures for T c = 0, 
77 = 0, and l/u = 0.39. 

To conclude, evidence of a zero-temperature transi- 
tion in the two-dimensional gauge glass with periodic 
boundary conditions is presented by studying the finite- 
size scaling of the correlation lengths. This approach is 
better than previous methods used because the correla- 
tion lengths show small statistical fluctuations and have 
a simple finite-size scaling form. In addition, the curva- 
ture of the spin-glass susceptibility as a function of sys- 
tem size indicates that there is no spin-glass transition at 
T « 0.20, at least for the system sizes studied. Moreover 
the curvature extrapolates within error bars to T c w 0. 
These results are in agreement with experimental data 
on YBCO thin filmsji^ where a finite-T is absent, and 
different numerical approaches In particular, it is 
shown that the system does not exhibit a spin-glass tran- 
sition at T « 0.20, as claimed in Refs. Hi, andE It is 
noteworthy that Refs.|8| and find the aforementioned 
finite temperature transition by studying the gauge glass 
with fluctuating twist boundary conditions instead of the 
commonly used periodic boundary conditions. Could this 
be one source of the discrepancies? 

Finite-size effects can influence the data, and so an 
analysis with larger system sizes is desirable. Nev- 
ertheless, the data presented here strongly support 
the absence of a transition at T « 0.20, although a 
transition at T < 0.13, the lowest temperature studied, 
while unlikely the case, cannot be ruled out completely. 
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